Columnar defects and vortex fluctuations in layered superconductors 
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We investigate fluctuations of Josephson-coupled pancake vortices in layered superconductors in 
the presence of columnar defects. We study the thermodynamics of a single pancake stack pinned 
by columnar defects and obtain the temperature dependence of localization length, pinning energy 
and critical current. We study the creep regime and compute the crossover current between line-like 
creep and pancake-like creep motion. We find that columnar defects effectively increase interlayer 
Josephson coupling by suppressing thermal fluctuations of pancakes. This leads to an upward shift 
in the decoupling line most pronounced around the matching fleld. 
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The question of the effective dimensionality of the vor- 
tex structure in highly anisotropic high-Tc superconduc- 
tors is the subject of intense experimental and theoret- 
ical interest. It was predicted that thermal fluctuations 
of vortex positions produce strong relative phase fluctu- 
ations in adjacent layers and that the phase coherence 
along the c-axis is destroyed above the decoupling line 
(DL) in the H~T plane Above this line an additional 
transverse field is not screened by Josephson currents and 
completely penetrates the sample. Several experiments 
were interpreted as a manifestation of decoupling ef- 
fects. 

In this Letter we focus on the influence of columnar 
defects produced by heavy ion irradiation on the vor- 
tex decoupling and the resulting effective dimensionality 
of the vortex structure. The interest in this question is 
motivated not only by the technological importance of 
highly anisotropic compounds such as Bi2Sr2CaCu20x 
(BSCCO), but also by the expected rich variety of non- 
trivial dynamic and static properties arising from the in- 
terplay between the continuous linear geometry of heavy 
ion tracks and the discrete pancake-like character of the 
vortex structure. 

An important question is whether columnar defects 
can be especially efficient to pin 2D pancakes which are 
nearly independent in each layers. Early experiments 
1^ of irradiation of BSCCO failed to show unidirectional 
pinning along the columns thus suggesting that the linear 
nature of columnar defects may not be efficient for highly 
anisotropic materials. However subsequent experiments 
revealed [||-|| that the relative enhancement of critical 
current (by several orders of magnitude) and the upward 
shift of the technologically useful irreversibility region, 
was even greater for highly anisotropic material such as 
BSCCO than for the less anisotropic YBCO. 

We demonstrate that the observed enhanced pinning 



of pancakes can be explained in terms of the recoupling of 
the layers by the columnar defects. We show that colum- 
nar defects effectively increase the interlayer Josephson 
coupling by suppressing thermal fluctuations of pancakes 
positions. This leads to an upward shift in the decoupling 
line, the maximal shift being expected around the match- 
ing fleld. We also investigate in detail pinning regimes 
of a single pancake stack on a single column and on the 
system of columns. 

The theory of vortex pinning by columnar defects in 
anisotropic 3D superconductors was developed in 
Mapping vortex lines onto the imaginary-time quantum 
mechanics of 2D bosons subject to static disorder leads 
to the prediction of a novel low temperature " Bose glass" 
phase with vortices localized near columnar pins, an in- 
finite tilt modulus and zero linear resistivity. The Bose 
glass transition line Bbg{T) is shifted to higher temper- 
atures as compared to the melting line B„i{T) of the 
clean lattice. The existence of the Bose glass phase 
is supported by several experiments A new fea- 

ture of highly anisotropic superconductors is the discrete 
nature of vortices. The physics of the system is gov- 
erned by three characteristic lengths: the vortex spac- 
ing a = (^o/B)^^^, the distance between defects d and 
the Josephson length 7s. Here B is the magnetic in- 
duction, $0 is the flux quantum, s the interlayer spac- 
ing and 7 is the anisotropy parameter. The Josephson 
length defines the crossover field Bcr = ^o/{lsY above 
which the interplane interaction energy of pancakes ex- 
ceeds the interlayer coupling. The discreteness changes 
the depinning temperature, localization length, binding 
energy, and critical current. We calculate these parame- 
ters using a discrete-time generalization of the mapping 
onto quantum mechanics. We show that at low temper- 
ature and moderate fields the system crosses over from 
pancake to line-like behaviour. We study this crossover 
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at short scales in thermodynamic quantities and in trans- 
port properties as a function of the current. At elevated 
temperatures and fields the layered nature of the super- 
conductor is essential and leads to a new phase, the "de- 
coupled Bose glass", where the phase coherence along 
the c-axis is lost but pancakes are aligned by columns 
and therefore the system exhibits infinite tilt stiffness. 

The statistical mechanics of a single stack of pancakes 
near a columnar defect is determined by the partition 
function: 
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Zjv(ro,rjv) / J]^ 

Q<n<N 
n— 1 



depending on the positions of the start and end points 
ro,r7v |l3j. The Josephson coupling binds the pan- 
cakes in the neighboring layers by a spring-like interac- 
tion with the elastic constant ei — (eo/7^) lii(7s/rcMt) 
where ei is the elastic tension in the continuous limit, 
eo = *&o/(4''''^)^i ^ is the London penetration depth for 
the currents flowing along ab planes and Vcut being the 
cut off radius of order of 5o if r„, r„_i < 26o and of order 
of |r„ — r„_i| otherwise. We assume 7 < A/s, and there- 
fore neglect the magnetic coupling as compared to the 
Josephson one. We model columnar defects as cylinders 
of normal material of the radius 60 • The pancakes inter- 
action with the columnar defect is given by an attractive 
pinning potential V{r), which at r — 60 S> ^ can be found 
from the London equations ||l0|] , 



F(r) = seo In (1-62/^2) 



(2) 



The maximum pinning potential per unit length is C/q ~ 
eo (ln(&o/C)+0.38). 

Using the standard transfer matrix method one finds 
that in the large- limit Zm{v,y') takes the form: 

Z^(r,r') = eM-E^sN/T -V{v')/T)c^o{v)<P^{v') (3) 

where Eq (the binding energy per unit length) and 0o (r) 
are the ground-state energy and the ground state wave 
function of the discrete analogue of a 2D "Shrodinger 
equation" : 



exp(-£;os/T)(/)o(r) 
eidv' 
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Taking the limit s — s- recovers the usual quantum me- 
chanical representation of vortex line. The probability 
of finding a pancake at point r near columnar defect is: 
P{v)^eM-V{r)/T)4>l{v) 

Localization effects are governed by the effective po- 
tential C/e//(r) = -(T/s) (exp[-T/(r)/r] - 1) which is 



characterized by the large dip in the close vicinity of the 
column and has a l/r^ tail far from the defect. The 
small distance dip dominates at low temperatures while 
the large distance tail becomes important at high tem- 
peratures (the exact criterion will be specified below) . At 
low temperatures the potential is effectively short-ranged 
and can be approximated by a delta- function Uef / (r) ~ 
-W6{v) with W = (27rT/s) ^.^^^ rdr(exp[-l^(r)/r] - 
1). The main contribution to W at low temperatures 
comes from short distances and can be estimated as 
{2-KThQrQ/ s)eyi\){sUQ/T), with the length tq being de- 
pendent on the detailed mechanism of the pancake nu- 
cleation at the column boundary. Approximating the 
pancake potential at distances < r — 5o ^ by 
V[r) w sC/o -f (seo/2) ln(l + {r ~ 5o) VC^) we obtain 
To ~ i{T / seoY^"^ . Within the delta-function approxima- 
tion Eq. (H) can be solved exactly giving the ground 
state energy as: 



So - (T/.s)ln{l -exp(-2^rV(eiM/))} 



(5) 



The distribution function is then P(r) = ?!'o('') + 
(sM^/T) (/)g(0)^(r), i.e., it consists of the delta-peak cen- 
tered at the defect having the weight 1 — exp(si?o/T) and 
the smooth part with the weight exp(si?o/r). A similar 
structure for P(r) has been observed in numerical sim- 
ulations . We define the depinning temperature T* 
by the condition that at T > T* the probability to find 
a pancake outside the column becomes higher than the 
probability to find it inside the defect. Matching proba- 
bilities of finding pancake inside- and outside the defect 
gives rise to the self-consistent equation 



T* = sC/o/ In 



sT* 
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where the quantities Uq and ei are taken at T = T*. 
Eq. ( H) is valid provided $;(T*) » bo, this is equiv- 
alent to the condition 7s ^ 60 and defines a range of 
parameters corresponding to the "pancake" regime of 
depinning as opposed to the "elastic string" regime de- 
scribed in Ref. Above T* the localization param- 
eters obey the "exponent-in-exponent" temperature de- 
pendencies, in particular, for the localization length 
= J (Pv r2(^g(r), one has 
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Evaluating the parameter W with the potential from (g) 
one finds that the l/r^ tail gives rise to a logarithmic 
divergence at large distances which has to be cut off 
at Finally one arrives at w ^T6ofoexp(^) + 

Tr&nCo In ^' \'i , where the first term is due to the short 
distances contribution. Substituting the expression for 
^; from Eq. (^ into the last equation we can see that 
the tail region dominates above the typical temperature 
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T** = 2sUo/\n{s'^Uo/blei) > T*. Numerical estimate 
for BSCCO with Aq = I8OOI, 7 = 200 60 = lOOl gives 
T* « 69 K, T** « 73 K. 

At temperatures significantly larger than T** discrete- 
ness is not important any more, Eq.( |^) reduces to the 
Shrodinger equation, and the localization parameters co- 
incides with those for the elastic string [TT|Jl^ 

= jsbo exp (^] ; Eo = exp (-7^] • (8) 

with Tstr — V^^^i^abo / TT . The single-column regim breaks 
down at temperature Tdi where the localization length 
matches the spacing between the columns, ^{Tdi) — d. 

Next we discuss transport of the pancake vortices in 
the presence of columnar defects at low temperatures and 
fields {B < Bcj, — ^o/d^ < Bcr) below critical current. 
The pancake motion in this region is controlled by the 
liberation of a single stack of pancakes from the localized 
state near a columnar defect. The energy cost to liberate 
a segment of N pancakes is: 

N 

AEn = (|^(r„ - r„+i)' + -^Uo + y(r„)) - ^jsr„ 

n—l 

The barrier V is determined by the lowest energy saddle 
point configuration determined by the stationarity con- 
dition SAEn /Svn = 0. One can distinguish the string- 
like and single pancake regimes separated by the "decou- 
pling" current jdcc = c^eiUo/ {^qs). At j < jdcc the 
critical nucleus consists of > 1 pancakes . The con- 
tinuum limit corresponds to A^ S> 1, i-e j jdec and 
one finds N ~ 2V2jdcc/j >> 1- The resulting string-like 
creep is characterized by a half loop width w ~ cUo/^oj 
and the barrier Vstr{j) ~ 4-\/2sC/ojdcc/(3j) and holds as 
long as w < a ,d. Because of discreteness the number 
of the pancakes in the nucleus decreases by one with in- 
crease of the current at well defi ned transition c urrents 

J = JN^N-I, JN^N-I = Jdcc\/2V2/N{N +1). Thus 

for jN+i^N < j < Jn^n-i, the barrier V{j) is deter- 
mined by the lowest energy for a N pancake excitation: 

Vij) = EnU) = Uo[N- ^N{N + l)(iV + 2)) . For 
j ^ jdec the last expression reduces to T^tr- Above jdoc 
pancakes liberate from the column independently with 
activation energy « sUq. 

We turn now to the collective localization of a pancake 
stack on the forest of columnar defects for B < B^ < 
B^ . There is a new "tight binding" regime for a localiza- 
tion of a single pancake stack by many columnar defects, 
which does not exist for continuous lines. Above the tem- 
perature Tj « tid^ / s pancakes can travel from column to 
column, but as long as T < T* one can neglect the prob- 
ability to find pancakes outside the columns. There are 
two sources for weak localization on the forest of identi- 
cal columns: fluctuation of column density and potential 



fluctuations due to superposition of individual long range 
column potentials. Here we focus on the density fluctua- 
tions which dominates at high temperatures, T > seobo/d 
and consider potential fluctuations elsewhere 0. The 
localization length 3> d can be estimated as follows. 
When confining the pancake stack within the region of 
size with enhanced number of columns N = N + 6N, 
N « {£,i/dY and 5N ^ \/N , we gain the positional en- 
tropy due to the enhanced number of available states as 
-{T/s){5N/N) w -{T/s){d/^i) . The entropy loss due 
to suppression of pancake thermal wandering by colum- 
nar defects is T^/ei^f. Optimizing the change in free 
energy 6F = -~{T/ s){d/^i) + T^/ei^f we find the local- 
ization length and binding energy: 

6«Ts/eid; Eb^-eiid/sf. (9) 

It is important to note that in this regime the effec- 
tiveness of localization decreases with increase of column 
concentration due to the decrease of elemental jump dis- 
tance. 

The above results allow us to estimate the typical cur- 
rent, jdoc ~ cEb/{^o£,i) ~ cefd^ / {^qs^T) at which the 
stack of pancake delocalizes from the weakly localized 
state. Note however that at T < T* the barriers control- 
ling vortex hopping between columns of the are still large 
(of order of s/7o), and thus for j > jdec the resistivity still 
remains much smaller than the flux flow resistivity p//: 
p « yO// exp(— ^^). The resistivity approaches p// at 
much larger current jc = cJ7o/($oO- The current jdec 
plays a role of "decoupling current" because at j > jdec 
the interlayer coupling does not influence much the mo- 
tion of pancakes. 

At elevated temperature pancakes spend more time 
outside the columns and at T = T* there is a crossover 
from the "tight binding" regime to the gaussian ran- 
dom potential localization. In the latter regime vortices 
are localized by the gaussian fluctuations of the effective 
potential (SC/e// « -(T/s) [exp (siJo/T) - 1] 60^0/(6^)- 
Adding the confinement term T'^/ei^f and optimizing 
leads to: 

C. = 4^(e^-ir;i^.--^(e^-l)^ (10) 
eiboro sd^ 

Eqs. (||,^0|) are valid at very small fields when inter- 
vortex interactions can be neglected. Interactions be- 
come relevant for B > Bcp ~ ^oEi,/{eo^f) and lead to 
a collective pinning regime . For the "tight binding" 
regime the crossover to collective pinning is expected at 

Bcp ~ BcrTj/T^. 

We now discuss the suppression of thermal decoupling 
by columnar defects. The effects of columns on the phase 
coherence between the layers is of a very different nature 
at high and low fields and have to be discussed sepa- 
rately. In the absence of columnar defects and at small 
fields B < Bcr = ^o/iis)"^ the decouphng occurs in the 
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liquid phase when the mean square fluctuation, uf^^^^ = 

^(rrt+i — of the separations between the adjacent 

pancakes belonging to the same vortex hne becomes of 
order 0. For pin- free system u^gji, = 2sT/ei and 
Tdeco(S) ~ ei$o/2sB ~ T™(Be./B)i/2 > Tra where r,„ 
is the melting temperature of the pure system. Colum- 
nar defects suppress the tilt deformations of vortex lines 
and reduce their influence on the interlayer Josephson 
coupling. Since the vortex-vortex interactions have little 
effect on u^^-^^ the above result for a single stack can be 
used. Correction to u^^jj^ caused by columns is connected 
to binding energy by relation Au^^jj^ = 2s'^{dEb/ dei). 
For temperature range T < T* decoupling takes place 
within the "tight binding" regime. Using we obtain 
2sT/ei — Const cP and the Lindemann-like crite- 



^ncib ~ ^"-^ /ti 

,2 ^ „2 



rion u^^i^ « a" gives 

Tdcc ~ Tdcco + Const ei(f/s, at Tdcc < T* 



(11) 



Note that in this regime the decoupling temperature 
decreases with the increase of the concentration of the 
columns due to decrease of the elemental jump distance. 
At Tdoc > T* one finds, similarly: 



dec ' 



dccO 



'"O'O 
d2 



exp 



Id 



- 1 



ccO 



(12) 



As temperature approaches T^, the renormalized decou- 
pling line smoothly merges the unrenormalized one. 

At large fields B 3> -B$ , -Bcr one has to take into ac- 
count that (i) only a small fraction of vortices w B/Bg, is 
trapped by columns, and (ii) the interaction between the 
vortex lattices in different layers is weaker than intralayer 
interactions. To estimate fluctuations of the phase differ- 
ence between layers we note that the phase disturbance is 
(5(/)„(q) — 2'!TnyUnt/<}±, where u„f is the transverse com- 
ponent of the lattice displacement in the n-th layer at the 
wave vector q. The resulting average phase difference be- 
tween the layers ((At/))^) = ((^0„(r) — (50„+i(r))2 
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((A0)2) = {2nn,f 



dq± {\unt{qi_W 
(2^)2- 



(13) 



The amplitude of the transverse vortex displacement can 
be estimated as (|u„t(q±)P) « T / {sC6e{q\+ {tt / df j) , 
where Cge = Beo/A^q is the shear modulus. Perform- 
ing integration we flnd ((A^)^) TB / {seoB^)ln{js / d) . 
Thus, the typical temperature Jdcc, above which the 
Josephson coupling is considerably suppressed can be es- 
timated as 



Tdc 



B^seo 
B\n (7s/rf)' 



(14) 



The ratio of Tdoc to the decoupling temperature for the 
ideal crystal Tdcco §, ccO 



($o/4^S)^/'eo/7, is 



Thus, at flelds B < B\/Bcr the columnar defects ef- 
fectively increase the coupled region. For B > B\/Bcr 
the curve Tdec smoothly merges with the decoupling line 
of the pure system. An important observation is that at 
large fields Tdec(-B) < and approaches T,„ at i? « 
The vortex phase in the region rdec(-B) < T < T„i has 
infinite tilt modulus and therefore can be referred to as 
the "decoupled Bose glass". The largest upward shift 
of the decoupling temperature is expected to occur at 
B « _B$ where the thermal motion of pancakes is maxi- 
mally suppressed and pancakes do not travel from column 
to column. 

To conclude, we have found three regimes of pan- 
cake localization on one column: low temperature, "pan- 
cake depinning", and "string depinning" regimes. A new 
"tight binding" regime of localization of a single pancake 
stack on many columns was found at high concentration 
of columns. We have studied the dynamic crossovers from 
line to pancake creep, which can be observed in transport 
experiments. We have shown that columnar defects "re- 
couple" Josephson layers by suppression of pancake posi- 
tional fiuctuations. This effect can be checked in neutron 
diffraction and /iSR experiments. We predict an upward 
shift of the decoupling line by columnar defects. A "de- 
coupled Bose glass" phase is predicted at high field. 
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